Domain wall effects in ferromagnet-superconductor structures. 



in 

O 

o 



Igor S. Burmistrov 1 ' 2 and Nikolai M. Chtchelkatchev 1,3 ' 4 'E 

1 L.D. Landau Institute for Theoretical Physics, Russian Academy of Sciences, 117940 Moscow, Russia 
Institute for Theoretical Physics, University of Amsterdam, 
Valckenierstraat 65, 1018XE Amsterdam, The Netherlands 
5 Institute for High Pressure Physics, Russian Academy of Sciences, Troitsk 14-2092, Moscow Region, Russia 
^Moscow Institute of Physics and Technology, Moscow 141 700, Russia 

We investigate how domain structure of the ferromagnet in superconductor-ferromagnet het- 
erostructures may change their transport properties. We calculate the distribution of current in the 
superconductor induced by magnetic field of Bloch domain walls, find the "lower critical" magneti- 
zation of the ferromagnet that provides vortices in the superconductor. 
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Superconductivity and ferromagnetism are two com- 
peting phenomena: while the first prefers antiparallcl 
spin orientation of electrons in Cooper pairs, the sec- 
ond forces the spins to be aligned in parallel. Their 
coexistence in one and same material or their interac- 
tion in spatially separated materials leads to a number 
of new interesting phenomena, for example, 7r-state of 
superconductor (S) - ferromagnet (F) - superconductor 
(SFS) Josephson junctions) 1 ' 2 ' 3 ! 4 ! 5 ! 6 ! 7 ! 8 highly nonmono- 
tonic dependence of the critical temperature T c of a SF 
bilayer as a function of the ferromagnet thickness 9 - and 
so on. Recent investigations of SF bilayers showed that 
their transport properties often strongly depend on the 
interplay between magnetic structure of the ferromag- 
net and superconductivityi 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 In par- 
ticular, it was argued that due to ferromagnetic domains 
vortices may appear in the superconducting film of the 
SF bilayer and the domain configuration, in turn, may 
depend on the vortices^ Recently, generation of vortices 
by magnetic texture of the ferromagnet in SF junctions 
was demonstrated experimentally^ In a number of ex- 
periments dealing with T c of SF bilayers, or Josephson 
effect in SFS structures 1 - the domain magnetizations were 
parallel to the SF interface. Ferromagnets used in the ex- 
periments were often dilute with the exchange field com- 
parable to the superconducting gap and small domain 
size [smaller or comparable to the bulk superconductor 
screening length] and broad domain walls. 1 

In this paper we try to put a step forward the answer 
to the question, how domain structure of the ferromag- 
net in SF heterostructures may change their transport 
properties. In the first part of the paper we discuss the 
junctions where S and F are weakly coupled (there is in- 
sulator layer in-between such that there is no proximity 
effect) and magnetizations of the domains are parallel 
to the SF interface. We find the distribution of current 
in the superconductor induced by magnetic field of the 
domain walls and the "lower critical" magnetization of 
the ferromagnet that provides vortices in the supercon- 
ductor. In the end of the paper we estimate the criti- 
cal temperature in strongly coupled SF bilayer when the 
proximity effect is strong. In this paper we do not con- 
sider the rearrangement of the domain configuration due 




FIG. 1: The SF junction. A sketch of a Bloch domain wall. 
The magnetization rotates according to Eq. Magnetiza- 
tion in the center of the domain wall is perpendicular to the 
S film. 



to the superconductor 1 ^ though we mention the super- 
conductor induced transitions between Bloch and Neel 
domain wall types. The point is that the crystal struc- 
ture of the ferromagnets in experiments 1 - was not perfect. 
Experimental data suggests that defects, dislocations in 
the lattice that appear during lithography process stick 
domain configuration. 

The domain texture of the F film is described by the 
following magnetization (see Fig. ^| 

oo 

M = MO (z) 8 (dp — z) (™l) J m(a;- jL). (1) 

j=— oo 

where vector m(x) rotates as follows 19 

7712 = 0, m y = tanh(x/<$), m z — — 1/ cosh.(a;/<5). (2) 

The vector potential A satisfies the Maxwell-Londons 
equation 

V x (V x A) + AX^O (-z) 6{z + d s ) = 4ttV x M (3) 

that should be supplemented by the standard bound- 
ary conditions of continuity of A and dA/dzm^ By 
solving Eq. J3J with a help of the Fourier transforma- 
tion we can find the distribution of the magnetic field 
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B = V x A in the entire spaced In agreement with 
general expectations the z component of the magnetiza- 
tion in the F film that collects at domain walls results 
in the current flow in the S film. It is convenient to de- 



fine the current averaged over the thickness of the S film, 
J y (x) — — c/(47rA|) f® d Aydz. Then we obtain 
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J y (x) = -2ncM 



sing„x 



<2 n sinh<3„ds + g„(cosh Q n d s - 1) 



*lL^ QnCOsh(lq n 5) (Ql 



q 2 .) sinh Q n ds + 2Q n q n cosh Q n ds 



[1 - exp(-q n d F )] 



(4) 



where c stands for the speed of light, q n = n(2n + 1)/L 
and Q n = \J q n + A^ 2 . Eq. (0J constitutes one of the 
principal results of the present paper. It allows to com- 
pute the distribution of the current flow in the S film for 
a general set of parameters ds, d F , S, L and A^. Below 
we shall analyze two the most interesting cases of thick 
(ds,dF ^> Xl) and thin (ds,dF <C Xl) SF bilayer. 

Thick SF bilayer. Eq. @ can be drastically simplified 
provided ds,dF ^> Xl such that the current J y {x) be- 
comes independent on the widths d$ and dp of the S and 
F films. It is given as 



J y (x) = -cM 



s—^ smq n x 
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2tt<5 

X\L ^ cosh %q n 6 Q n (Qn + q,i 



(5) 



In order to understand the distribution of the current 
J y (x) as determined by Eq. (J5J we shall first analyze the 
case of a single domain wall. Taking the limit L — > oo in 
Eq. JSJ, we obtain the following result for the current in 
the presence of a single domain wall in the F film 



Jy( x ) 
cM 



S 



du> 



1 



j 2 cosh 



(6) 



The distribution of the current J y (x) is governed by the 
single parameter it6/(2Xl) as it is shown in Fig. [21 If 
the width S of the domain wall is small compared to the 
Londons penetration length A^, 7t<5/(2Al) -C 1, we find 
the distribution of the current as 



J y (x) = cM 



- 1 
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\x\ > 6, 
\x\ < 5, 



(7) 



where Ki(x) is the modified Bessel function of the sec- 
ond kind. In the opposite case of the thick domain wall, 
7t<5/(2Al) ^> 1, we obtain 



J v (x) , x 2 , / 1 . x 
yy ' = tanh Im ib - + i — - 



Xl tanh | 
(5 cosh | 



(8) 



where VK 2 -) denotes the digamma function. 

According to Eqs. and JHJ) the current J y (x) be- 
haves linearly with x for x *C <5 and decays as a power 
law for large x. The current distribution J y (x) is spread 



on the length L s oc max{<5, Al} from the origin whereas 
its maximal value J™ oc cMS/L s = cM min{ 1, S/Xl}- 

Now we turn back to the general case of multi domain 
wall structure in the F film that corresponds to a finite 
size L of domains. We have evaluated the sum in Eq. (J3J) 
numerically and present results for the current distribu- 
tion in Fig. |3 While Xl remains small compared with 
L the profile of J y (x) corresponds to almost independent 
current distributions near each domain wall that results 
in distinctive two maximum structure as shown in Fig. [3] 
When Xl becomes of the order of L the two maximum 
structure transforms into sinusoidal profile with the max- 
imum exactly in the middle of a domain. 

Thin SF bilayer. In the case of the thin SF bilayer, 
ds,d,F -C Xl, by expanding the general expression (@J in 
powers of ds and d F , we find 



J y (x) 



-cMdt 



2tt5 



srng„x 



coshf q n S 1 

n— z 



2q n X 



(9) 



Here A = X\/ds usually referred as the effective penetra- 
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tion lengthi*! We shall first analyze the case of a single 
domain wall again. In the limit L — > oo we obtain from 
Eq. © 



, . d F S 
J y (x) = -cM^ 



duj- 



u> cosh 2io 



(10) 
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The distribution of the current J y {x) is presented in 
Fig. 01 for different values of the parameter tt5/(4X). 
If the domain wall is thin, 7n5/(4A) < 1, Eq. (|TJ3l yields 



Jy( x ) 
cM 

where 
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|x| > 6, 
|x| < (5. 



(11) 



(12) 



Here /(x) = sinxci(x) — cosxsi(x) with ci(x) and si(x) 
being the cosine and sine integral functions. In the op- 
posite case 7T(5/(4A) 3> 1 the current distribution J y (x) is 



given as 



J y {x) dp tanhx/<5 
cM 5 coshx/<5 



(13) 
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FIG. 2: The case d > X L . The plot of J y (x)/(cM) as a 
function of x/X L for ir6/{2X L ) = 0.1, 0.5, 1, 2, 4. 



FIG. 4: The case d < A L . The plot of 2XJ y (x)/(cMd F ) as a 
function of z/(2A) for tt5/(4A) = 0.2, 0.5, 1, 2, 4. 
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FIG. 3: The case d > X L . The plot of J y (x)/(cM) 
as a function of x/L for 5/Z = 0.02 and Xl/L — 
0.005,0.01,0.05,0.1,0.25,0.5. 



Eqs. (|TTT) and (|T3)l proves that the J y (x) increases 
linearly with x for x <C 5 and decreases algebraically 
for large x. The position of the maximum of J y (x) 
is situated at L s oc £ and the value at the maximum 
Jy l oc (cLf/6) min{l, S/X}. As one can see therefore the 
current distribution for the thin SF bilayer is qualita- 
tively different from one for the thick SF bilayer. 

In the case of multi domain wall structure in the F 
layer we have performed evaluation of the sum in Eq. @ 
numerically and have obtained the results for the current 
distribution J y (x) presented in Fig. [SJ We mention that 
the two maximum structure survives even for A of the 
order of L for the thin SF bilayer. 

As known the lower critical field for a thin S film is 
much smaller than for the bulk superconductor. There- 
fore it is possible that even small magnetization collects 
at domain walls can induce a vortex in the thin S film. 12 
Let us assume that there is a single vortex in the the thin 
S film situated at x — 0. The magnetic field becomes a 
sum of the magnetic field induced by the domain walls 
and the magnetic field of the vortex. The free energy can 
be written as 

i^ = yd 3 r^ + ^|VxB| 2 -MBj, (14) 
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FIG. 5: The case d < X L . The plot of J y (x)/(cM) as a 
function of x/L for S/L = 0.02, d F /L = 0.1 and X/L = 
0.05,0.1,0.25,0.5,0.75, 1. 



where B denotes the total magnetic field. The difference 
T of the free energy for the state with the vortex and 
the free energy for the state without vortex is given as 
follows 2 ^ 
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' Mc = ^~d7 g{6 ' x > L) > (15) 



where 0o = ch/(2e), H Cl = (0 o /47rA 2 ) In A/£ is the lower 
critical field in the thin S film without the F film and 
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i 



2tt 



L ^ cosh ^ (l + 2g„A) 2 



(16) 



The T becomes negative if M > M c and vortices can pro- 
liferate in the S film until vortex- vortex interaction stops 
it or (that is more probable) the domain wall changes to 
Neel domain wall type to reduce the free energy. In the 
most interesting case of a single domain wall we find 



7r<5/4A < 1, 



(17) 



4tt d F 1 l-32GA/(7r 2 <5), nS/4X » 1 
where G w 0.916 stands for the Catalan constant 
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(Note, that M c given by Ea. H17(l differs from estimates 
of Ref[H). 

In conclusion, domain wall effects in ferromaget- 
superconductor structures are investigated. We find the 
distribution of current in the superconductor induced by 
magnetic field of Bloch domain walls, calculate the "lower 
critical" magnetization of the ferromagnet that provides 
vortices in the superconductor. 

We neglected above the proximity effect in SF struc- 
ture assuming that the ferromagnet and the superconduc- 
tor are weakly coupled. Below we discuss the case when S 
and F are strongly coupled. Consider a SF bilayer with 
a perfect SF boundary. When the superconductor and 
the ferromagnet are thin enough then the bilayer can be 
described as a "ferromagnetic superconductor" with ef- 
fective parameters^ the superconducting gap A e g, the 
exchange field Ec^... The superconductivity survives 



in this system until E C x < A^ ff , where A^ ff ' is the gap 

at -Eex^ = 0. 22 Domain wall structure of the ferromag- 
net makes the effective exchange field nonhomogeneous. 
We find that if E^' changes its sign on scales of the 
order of £o or smaller then superconductivity in the bi- 
layer survives at yj ((E%? ] ) 2 ) > A<$, where ((E& S) ) 2 ) 
is the average square of the effective exchange field over 
the sampled 
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